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ABSTRACT 


The  basic  statistical  properties  of  the  wall  pressure 

field  associated  with  turbulent  spots  in  a transition  boundary 

layer  have  been  measured  on  a large  flat  plate  test  fixture 

in  the  Anechoic  Plow  Facility  at  the  Naval  Ship  Research  and 

Development  Center.  The  measured  statistical  properties  and 

the  distribution  of  the  turbulent  spots  on  the  plate  surface 

are  compared  with  a quantitative  formulation  of  the  transition 

1 2 

process  which  was  developed  by  Emmons  ' . An  important 
unknown  parameter  in  Emmons'  model  of  the  transition  process 
is  the  burst  source- rate  density  function.  This  paper  con- 
siders the  consistency  between  the  predicted  and  measured 
statistical  burst  properties  when  two  simple  forms  of  the 
source-rate  density  function  are  assumed.  Experimental 
values  of  the  burst  rate  and  average  burst  length  are  presented 
as  a function  of  the  intermittency  factor.  The  statistical 
results  are  normalized  in  term  of  spatial  and  temporal 
individual  burst  properties  and  compared  with  the  predicted 
results. 
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INTRODUCTION 


When  a streamlined  body  moves  at  high  speed  through  a viscous 

fluid  the  laminar  boundary  layer  flow  along  the  surface  of  the 

body  usually  undergoes  transition  to  turublent  flow  at  some  location 

on  the  forward  nose  section  of  tne  body.  Between  the  laminar  and 

fully-developed  turbulent-flow  regions  of  the  boundary  layer,  a 

transition  region  exists  over  a finite  length  of  the  nose  i which 

the  flow  is  observed  to  be  intermittent,  that  is,  patches  of 

turbulent  flow  occur^  randomly  within  the  laminar  boundary  layer 

3 

region  on  the  surface.  Detailed  measurements  have  shown  that  these 
turbulent  bursts  grow  as  they  are  swept  downstream  until  they  merge 
together  forming  the  fully-developed  turbulent  boundary  layer. 

An  understanding  cf  the  parameters  that  govern  the  width  of 
this  transition  region  along  with  detailed  data  on  the  convection 
and  growth  properties  of  the  turbulent  bursts  is  essential  to  the 
formulation  of  an  appropriate  wall  pressure  field  forcing  function 
required  for  the  prediction  of  the  flow  induced  noise  which  is 
transmitted  through  the  structures  of  the  moving  body. 

The  basic  properties  of  the  wall  pressure  field  associated 
with  turbulent  spots  have  been  measured  in  detail  on  a large  flat 
plate  in  the  Anechoic  Flow  Facility^’*\'AFF)  at  the  Naval  Ship 
Research  and  Development  Center  (NSRDC). 
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This  paper  concentrates  on  describing  the  statistical  properties 
and  distribution  of  turbulent  spots  during  natural  transition  on  a flat 
plate  and  presents  a statistical  model  of  a source  rate  density  function 
initially  developed  by  Emmons^. 


EMMONS'  MODEL 

In  1951, Emmons  first  advanced  the  concept  of  turbulent  burst 
generation  in  a laminar  boundary  layer  from  his  observations  in  a 
water-table  flow  experiment.  From  visual  observations,  he  concluded 
that  the  randomly  generated  turbulent  bursts  act  independently  of  each 
other  and  uniformly  grow  as  they  are  swept  downstream  by  the  flow. 

In  order  to  develop  a quantitative  formulation  of  the  transition  pro- 
cess, Emmons  assumed  the  existence  of  a source-rate  density  function. 
This  function,  g(x0,yQ.t0),  specifies  the  rate  of  production  of 
turbulent  point-source  bursts  per  unit  area  on  the  surface  and  is 
dependent  on  the  position  on  the  surface  (xQ,y0)  and  time  of  occur- 
rence (tQ). 

Emmons  developed  a probabilistic  model  which  related  the 
statistical  properties  of  turbulent  bursts  for  any  prescribed  form 
of  g(*0 «VV*  deve^°Pment  °f  Emmons'  model  was  based  on  prob- 
ability theory  and  his  limited  data  on  the  growth  characteristics  of 
turbulent  bursts.  New  data,  obtained  in  Reference  (3),  have  sub- 
sequently clarified  the  growth  properties  since  Emmons'  observations 
with  the  simple  water-table  experiment.  A brief  outline  of  the  con- 
cepts involved  with  the  formulation  of  this  model,  as  presented  by 
Emmons,  will  now  be  given. 
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Consider  a point-source  p0(x0»y0»t0)  f°r  the  creation  of  a burst 
on  the  body  surface.  In  che  x,y,t  space,  this  burst  sweeps  out  a cone- 
like  volume  whose  shape  depends  upon  the  convective  and  growth  pro- 
perties of  individual  bursts.  This  volume  which  originates  upstream  at 
a point  source  PQ  is  called  the  Propagation  Cone  and  is  the  domain  of 
influence  of  a source  at  p0(*0»y0*t0). 

Now  consider  a point  P(x,y,t)  at  which  a Flush-mounted  pressure 
transducer  is  positioned  to  monitor  the  intermittent  wall-pressure  field. 
The  time  record  of  this  pressure  field  signature  will  indicate  intermit- 
tent patches  of  boundary  layer  turbulence.  The  fraction  of  the  total 
sample  time  that  the  boundary  layer  at  a given  point  P is  turbulent 
is  defined  as  the  intermittency  factor,  y.  In  the  computation  of  the 
intermittency  factor,  one  must  add  up  the  burst  periods  (xg)  from  all 
bursts  originating  upstream  that  pass  the  point  P in  the  sample  time 
T total*  However,  one  must  exclude  those  portions  of  overlapping  bursts 
that  originate  at  different  upstream  locations  but  coincide  within  the 
sampling  time. 

The  locus  of  all  points  PQ  which  can  influence  the  state  of 
turbulence  at  point  P is  a volume  upstream  called  the  Retrograde  Cone 
(R).  This  cone-like  volume  has  propagation  rays  drawn  upstream  from 
P reflecting  the  domain  of  dependence  of  point  P.  All  bursts  which 
pass  over  the  point  P must  have  originated  within  the  dependence 
volume  R.  Figure  1 shows  the  projection  of  the  Propagation  and 
Retrograde  Cones  on  the  plate  surface. 

From  the  basic  concepts  of  the  growth  and  convective  properties 
of  turbulent  bursts  it  is  possible  to  use  probability  theory  to 
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analytically  derive  certain  statistical  relationships  between  the 
transition  region  and  the  properties  of  the  turbulent  bursts.  Three 
such  statistical  relationships  analytically  derived  by  Emmons  are  the 
Intermlttency  factor,  y,  the  average  burst  period,  ig,  and  the  burst 
rate  frequency,  fg,  which  Is  the  frequency  at  which  bursts  pass  a 
given  point.  The  derivation  of  these  statistical  relationships  will 
now  be  given. 

INTERMITTENCY  FACTOR 

The  probability  that  a burst  originating  at  PQ  will  be  detected 
by  a downstream  transducer  located  at  P is  defined  by  Emmons^  as 


♦<P.P0)  = 

If  we  let 


1 if  PQ  is  located  in  Region  R of  P for  x>xQ 

0 if  P„  is  not  in  Region  R of  P for  x>x„ 
o o 


dA„  = dx  dy„  and  dA  = dxdy 

o o ’'O 

dV  = dx  dy  dt.  dV  = dxdydt, 

O O O 0 

then  we  would  be  tempted  to  say  that  the  fraction  of  tine  that  P is 
turbulent  is  given  as 


y(P)  - f <)>(P,P  )g  (P.)  dV  (1) 

JR 

where  the  region  of  integration  of  4>(P,P0)g(P0)dV0  is  over  the 
dependence  volume  R.  However,  note  as  previously  mentioned  that  more 
than  one  burst  can  be  generated  within  the  Region  R such  that  they 
will  coincide  in  time  t at  the  transducer  position  P.  For  example. 
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bursts  originating  at  P^O^.t*)  and  P0(x0»y0»t0)  "hen  x‘  < xQ  and 
t^  < can  be  detected  colncldentlally  (see  Figure  1).  For  this  case, 
two  such  bursts  coinciding  at  a point  P and  time  t would  be  counted 
separately  giving  twice  the  actual  value  of  intermlttency.  To  avoid 
these  burst  overlaps,  It  Is  assumed  that  the  burst  originating  nearest  the 
leading  edge  of  the  plate  will  be  counted  as  causing  the  turbulent 
burst  of  length  tg  at  the  transducer  position  P since  It  will  be  of 
longer  duration  and  have  larger  burst  length.  This  Is  true  provided 
that  the  burst  growth  rate.  Ug,  and  burst  convection  velocity,  Uc,  are 
Identical  for  both  bursts  alono  their  downstream  paths. 

Now,  let  ip(P ,F0)  dVQ  be  the  probability  that  a turbulent  burst  is 
detected  at  point  P due  to  sources  ii?  dVQ  at  a point  PQ  but  not 
because  of  any  sources  at  PQ(xQ»y‘.t*)  with  x*  < xQ,  since  the  source 
at  PQ  Is  to  be  counted  only  If  it  is  the  point-source  nearest  the  lead- 
ing edge.  From  this  probability  function,  Emmons  concludes  that 

*(P.P0)  = 0(p,p0)  g(P0)dV0][l  - f U>(P,Pj)dV»]  (2) 

R1 

where  g(P0)  is  again  the  source-rate  density  function  and  Rl  Is  the 
region  shown  in  figure  1.  The  first  term  in  Equation  (2), 
[<J>(P.P0)g(P0)dV0],  represents  the  fraction  of  time  that  bursts  will  be 
detected  at  P due  to  sources  at  PQ.  The  second  term  [1  - j iHP.P^dV*], 
is  the  fraction  of  time  that  no  bursts  are  observed  at  P dul  to  sources 
upstream  of  xQ.  Thus  the  intermittency,  or  fraction  of  the  total 
sample  time  that  bursts  are  being  detected  at  P is  given  by 
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^(P,P  ) <N 
o o 


(3) 


Y(P) 


R 


Emmons  solved  the  integral  equation  obtained  by  substituting 
Equation  (2)  into  (3)  and  found 


Y(P)  - 1 - exp[-j  g(PQ)  dvj  (4) 

R 

Equation  (4)  gives  the  relationship  between  the  intermittency  factor 
and  the  form  of  the  source-rate  density  function  g(P0) . This  relation- 
ship is  solely  dependent  on  the  form  of  g(PQ)  and  the  shape  of  the 
dependence  volume  R. 

BURST  RATE  FREQUENCY 

During  a sample  time  Tt  • point  P will  be  turbulent 

Ttotal  ’KP.Pq^o  seconds  because  of  sources  in  dVQ.  The  number  of 

bursts  passing  P from  sources  at  P„  in  dVA  is 

oo 

*<P-Po>  dVo  . W"''.1''. 

^ (5) 
max 

where  xg  is  taken  as  the  burst  length  measured  at  P(x,y,t)  due  to  a 
burst  originating  upstream  at  P (x  ,y  ,t)  and  xg  is  the  maximum 
burst  length  along  the  centerline  of  the  burst.  When  y = yQ,  the 
burst  erupts  along  an  upstream  streamline  and  the  sampled  burst 
length  xg  equals  x„  . However,  for  the  general  situation,  one  can- 
not  distinguish  between  xg,  the  burst  period  along  an  arbitrary 
section  of  the  wedge-shaped  burst  and  xg  . Since  xg  is  the  only 
quantity  amenable  to  an  analytical  formulation,  the  equations  are 
expressed  as  inequalities.  Reference  (4)  includes  data  on  the 
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distribution  of  burst  periods  obtained  from  observations  at  a fixed 
position.  The  total  number  of  bursts  per  second  passing  P that 
originate  In  the  domain  of  dependence  region  R Is  then 


N(P)  - T 


total 


f »<p-p„>  <*„  > f 
Jr  tb  ~ 'r  tb 

max 


(6) 


Thus  the  frequency  at  which  bursts  pass  point  P Is  the  burst  rate 
frequency  given  by 


urn  f ^P*PJ  dVft  f w,* 

fB(?)  - > - 

B Ttotal  J TB  Jr  t 


^(p,p0)  dVo 


(7) 


'R  B 


sax 


In  terms  of  the  source- rate  density  function  Equation  (7)  Is  express- 
ed In  the  form 


•1- 

Jr 


g(P  )exp[ 

fB(P)  - I 0 JR 


-[  g(P*)  dV> 

Jr*  0 0 


I 


dV 


‘B 


-JE 


g(P^)exp[-j  ^(P^dV*] 


dV 


(8) 


B 


max 


AVERAGE  BURST  PERIOD 

The  average  burst  period  (length)  Is  given  by 


Y(P)  Ttptal  y(p) 

• T1 


B N(P) 


fB(p) 


(9) 


and  can  be  determined  once  the  values  of  y and  fg  are  known. 

As  can  be  seen  from  these  equations,  the  statistical  relationships 
expressed  In  Equations  (4),  (8),  and  (9)  are  functions  of  only  the 
source  rate  density  function,  g(PQ ) , and  the  growth  and  convective 
properties  of  turbulent  bursts.  If  the  source- rate  density  function 
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is  known  and  the  growth  and  convective  burst  properties  are  measured, 
it  Is  possible  to  find  at  least  a numerical  solution  to  Equations  (4), 
(8),  and  (9). 

In  the  solution  to  Equations  (4),  (8),  and  (9)  we  may  assume,  for 
our  simple  case  of  a steady  mean  flow  over  a flat  plate,  that  the  flew 
properties  are  functions  of  only  the  streamwise  coordinate  x.  This 
assumption  gives 

8(p0)  - 8(*0.y0.te)  - g(xo)  (10) 

The  Integrations  In  Equations  (4)  and  (8)  *re  performed  over 
the  domain  of  dependence  R,  thus  the  Incremental  volume  element  dV 

o 

can  be  expressed  In  terms  of  the  growth  and  convective  properties  of 
the  individual  bursts.  This  transformation  of  the  incremental  volume 
element  Is  shown  In  Appendix  A and  is  given  as 


(x-x  ) 2 

dV  ■ — a*  dx 

o U o 


(ID 


and  xD  is  given  as 
Bmax 


T 


B 

max 


U 

U U 


(lib) 


where 

Ufi  = freestream  velocity  at  the  edge  of  the  boundary  layer 


u « 

°*  " uu  tan(0t)  (lie) 
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(lid) 


Ug  * VUt*  UC  S 1/2  {W 

Ut  3 trailing  edge  velocity  of  burst 
tan  ( «)  * half  angle  of  burst  (see  Figure  2) 
and  = leading  edge  velocity  of  bursts 

Equations  (4)  and  (8)  then  become,  respectively 

Y(x)  - 1 - exp  [—  f g(xQ)(x-xo)2  dxQ]  (12) 

e 'o 


and 

fB(x)  > tan  (a)  fx  (x-xo)g(xQ)  exph^p  [ °g(x*> (x-x*)2  dx‘]  dxQ  (13) 

To  apply  Equations  (12)  and  (13)  to  the  transition  on  a flat 
plate  the  form  of  the  source-rate  density  function  o must  be  assumed 
or  deduced. 

There  are  two  different  approaches  by  which  g(xQ)  in  Emmons' 
model  may  be  deduced  from  experimental  iata.  Thii  first  approach  is 
to  successively  differentiate  Equation  (12)  to  obtain  the  relationship 
given  by  Narasimha/, 


g(x) 


u 

e 

2o* 


d3G 

dx5 


(14) 


where  6 = -ln(l-y).  From  this  relationship  measured  values  of  the 
intermittency  distribution  can  be  used  to  determine  the  form  of  g(x). 

The  second  approach  is  to  assume  various  forms  of  g(x)  and  solve 
Equations  (12)  and  (13).  These  analytical  results  can  then  be  compared 
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to  experimentally  measured  values  of  y(x),  fg(x),  and  Tg(x).  This 
process  allows  the  determination  of  that  form  of  g(x)  which  most 
closely  predicts  the  turbulent  bursts  properties. 

It  must  be  realized  that  in  the  first  approach,  the  third 
derivative  of  an  experimentally  determined  distribution  must  be  taken. 
The  merits  of  such  an  operation,  in  view  of  the  accuracy  of  the  data, 
is  highly  questionable.  Thus  the  more  realistic  approach  Is  to  use  the 
assumed  forms  of  g(x)  to  solve  the  analytical  expressions  and  then 
compare  these  results  to  the  experimental  results. 

SOLUTIONS  FOR  PARTICULAR  SOURCE  FUNCTIONS 

Enmons1  inodel  will  now  be  solved  for  two  assumed  forms  of  g(x). 
CASE  1 - q(x)  DELTA  FUNCTION 

The  first  assumed  form  of  g(x)  is  that  of  a Dirac’s  Delta 
function. 


g(x)  - n<5  (x-xt)  (15) 

where  n is  defined  as  the  number  of  sources  per  unit  length  per  unit 
time  along  a line  located  at  x - xt  from  the  flat  plate's  leading 
edge  and  is  dependent  on  the  freestream  velocity  Ug.  The  physical 
meaning  of  this  form  of  g(x)  is  that  all  sources  of  turbulent  bursts 
are  located  along  a line  perpendicular  to  the  flow  stream.  Sub- 
stitution of  Equation  (15)  into  (12),  (13)  and  (9)  results  in  the 
relations 

Y - 1 - exp  b—  C ] 0^) 

l e 
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(17) 


fB  > «Bl  " ntan(a^t  exp[^“  £*] 


T . Y(x)  < - . 

TB  fB(x)  - TBj  fB  (x) 

l 


(18) 


where  Is  defined  as  x-x^. 

Solving  Equation  (16)  for  In  terms  of  Yt»  one  obtains 


(19) 


Equation  (19)  gives  a relationship  that  defines  the  distance  down- 
stream from  the  transition  point  at  which  a given  value  of  intermit- 
tency  occurs.  This  relationship  can  be  used  to  find  the  width  of  the 
transition  bursting  region,  AXt>  where  AXt  is  defined  as 


Ax  - (x-x  ) - (x-x  ) 

1 Y-0.99  Y-0.01 

By  substituting  Equation  C9)  into  (20)  we  get 


(Ax  ) 
t i 


(20) 


(21) 


where  it  is  assumed  that  Ug  is  essentially  constant  over  the  transition 
region. 

CASE  2 - g(x)  UNIFORM  SOURCE  DISTRIBUTION 

As  a second  simple  assumed  form  of  g,  g is  taken  as  independent 
of  position  on  the  plate.  This  condition  is  expressed  by 
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g(x)  » g ■ constant. 


(22) 


where  g Is  defined  as  the  number  of  sources  per-unit  surface  area  per 
unit  time.  Although  g Is  independent  of  x It  Is  dependent  on  the  free- 
stream  velocity  U . The  physical  meaning  of  this  form  of  g(x)  Is  that 
turbulent  burst  sources  are  uniformly  distributed  within  the  transition 
region.  Substituting  Equation  (22)  into  (12)  one  obtains  the  relations 


Y 

2 


1 - exp[ 


-go* 


f*  1 

3U  V 
e 


gtan  (a)  exp[ 


_ga* 

3U 


qi 


r* 

* n 


; exP (3U”  c’l 


(23) 

(24) 


T 


B 


ll*l  < 7 . 1 1*1 
fB(*>  ~ \ fB  (*) 

2 


Solving  Equation  (23)  for  we  find  that 


(25) 


x-x 


t 


(26) 


and  substituting  for  values  of  x at  which  y = .99  and  { = .01  we 
J 1 2 2 

find  that 


(Ax.) 
t 2 


1.448 


(27) 


where  Ufi  Is  again  assumed  to  be  constant  over  the  streamwlse  transition 
width. 
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In  order  to  predict  the  width  of  transition  using  either  Equation 
(21)  or  (27),  the  relationship  between  the  source  rate  density  g and 
the  free st re an  velocity  Ug  must  be  known. 

EXPERIMENTAL  RESULTS 

Measurements  were  made  of  the  spatial  and  temporal  properties  of 
turbulent  pressure  bursts  In  the  boundary  layer  transition  region  on 
a flat  plate  (Appendix  8 discusses  the  method  by  which  the  spatial 
and  temporal  properties  of  the  turbulent  pressure  bursts  were  measured). 
The  empirical  results  presented  i.i  this  section  will  be  compared  with 
those  predicted  by  Emmons'  model  in  the  following  section. 

Figure  2 shows  the  three-dimensional  shape  of  Individual 

3 

turbulent  bursts  determined  by  Schubauer  and  Klebanoff  using  hot  wire 
anemometers.  Ti.e  values  of  the  burst  parameters  shown  in  the  figure 
were  those  determined  in  the  present  study  using  flush  mounted  pressure 
transducers. 

Figure  3 shows  the  experimental  intermittency  distributions 
through  the  transition  region  in  terms  of  a normalized  Reynolds  number. 

Figure  4 presents  the  measured  values  of  burst  frequency,  fg, 
plotted  versus  y.  Faired  lines  are  drawn  through  those  data  points 
which  represent  measurements  made  at  the  same  downstream  coordinate. 

Figure  5 gives  the  measured  mean  burst  length  as  a function  of 
the  Intermittency  factor. 


14 


COMPARISON  OF  EXPERIMENTAL  RESULTS  WITH  EMMONS'  MODEL 


in  order  to  directly  compare  the  experimental  results  with  those 
predicted  using  the  different  forms  of  the  source  rate  density 
function,  g(x),  in  Emmons'  model,  the  measured  data  are  normalized 
relative  to  the  maximum  burst  frequency  (the  normalization  details 
are  given  m Appendix  C).  Figures  6 and  7 show  the  comparison  between 
the  experimental  data  and  the  theoretical  results  obtained  using, 
respectively,  a Dirac  Delta  function  and  a constant  for  the  source 
rate  density  function  in  Emmons'  model. 

As  shown  in  Appendix  C,  a functional  relationship  between  the 
freestream  velocity,  Ug,  and  the  magnitude  of  the  source  rate  density 
function  can  be  determined.  Figures  8 and  9 show  these  relationships 
for  the  two  cases  considered.  Also  shown  in  the  figures  are  values 
of  the  width  of  the  transition  region,  Ax^,  calculated  using  the  n 
versus  Ufi  (Equation  (Cl 3)) and  g versus  Ug  (Equation  (C23) ) relationships 
and  Equations  (21)  and  (27),  respectively.  The  five  values  of  Axt 
shown  on  each  figure  were  calculated  from  the  corresponding  data 
points  relating  g and  Ug  which  were  calculated  using  Equations  (C13) 
and  (C23)  of  Appendix  C. 


DISCUSSION  AND  CONCLUSIONS 

? 

Emmons  developed  a probabilistic  model  of  the  statistical 
burst  properties  for  the  general  case  of  a source-rate  density  g(x) 
being  dependent  on  the  downstream  spatial  location.  However,  lack 
of  available  experimental  data  did  not  permit  him  to  examine  carefully 
the  explicit  nature  of  this  g(x)  function  except  for  the  case  of  g 


equal  to  a constant.  Emmons  attempted  to  evaluate  the  constant  g by 
using  a few  1 /5-second  strips  of  hot-wire  records  published  by 

O 

Schubauer  and  Skramstad  . Similar  calculations  were  made  by  the 
authors  using  the  new  data  included  in  this  report.  The  results  are 
shown  in  Figure  9. 

Narasimka7  assumed  the  form  of  g(x,.y„,tj  to  be  a Dirac  Delta 
function  and  evaluated  F.mmons'  equations  for  the  resulting  y(x) 
distribution.  From  the  comparison  between  experimental  data  and  the 
theoretically  obtained  y(x)  distribution,  Narasimha  concluded  that  the 
turbulent  burst  generation  process  does  indeed  occur  like  a delta 
function  source  along  a line  transverse  to  the  flow.  It  has  been 
found  by  the  authors  that  the  y(x)  distribution  for  the  constant  source 
model  is  very  similar  to  the  y(x)  distribution  for  the  delta  source 
model.  The  data  shown  in  Figure  3 was  also  fitted  to  a Gaussian 
Integral  Curve.  The  lack  of  sensitivity  in  the  y(x)  distribution  to 
the  various  forms  chosen  for  g(x)  leads  one  to  conclude  that  the  y(x) 
distribution  should  not  be  used  as  a criterion  for  determining  the  g(x) 
function. 

It  should  be  noted  that  in  the  present  experiment  the  streamwise 
extent  of  the  transition  region  coui-  not  je  accurately  measured  at 
a given  flow  velocity  due  to  the  use  of  a limited  number  of  fixed, 
wall-mounted  pressure  transducers  for  detecting  the  wall  ^essure 
bursts  in  the  intermittent  transition  region.  Rather,  at  a given 
streamwise  position  the  ' ntermi ttency  was  varied  from  burst  onset 
through  fully  turbulent  fiow  by  varying  the  flow  velocity.  As  shown 
in  Figure  3,  the  use  of  the  normalized  Reynolds  number  results  in 


good  collapse  of  the  experimental  intermittency  distributions  determined 
through  the  transition  region  by  varying  flow  velocity  at  different 
transducer  locations  on  the  plate  surface. 

An  important  characteristic  of  the  transition  process  is  the 
point  at  which  the  burst  frequency  maximizes.  In  Figures  6 and  > this 
point  can  be  seen  to  be  at  intermittency  values  of  approximately  0.4 
and  0.8  for  the  Dirac  Delta  source  and  the  constant  source,  respectively. 
The  maximization  point  in  the  experimental  burst  frequency  curves 
occurs  at  an  intermitter.cy  value  of  approximately  0.6.  This  experi- 
mentally determined  maximum  point  falls  between  the  two  analytically 
determined  values  for  the  uniform  and  line  burst  source  distributions. 

As  seen  in  Figures  6 and  7 the  normalized  experimental  data  for 
the  case  when  the  burst  generation  occurs  at  a line  source  is  in  better 
general  agreement  with  Emmons'  theory  than  for  the  case  of  a uniform 
source  distribution. 

Inherent  in  any  analytical  model  for  the  burst  source  distribution 

is  a method  for  predicting  the  streamwise  width  of  the  transition 

region.  The  predicted  width  of  the  transition  region  given  by  Equations 

(21)  and  (27)  is  seen  to  depend  on  the  freestream  flow  velocity  Ug 

(Ufi  is  assumed  to  be  constant  over  the  transition  region  for  the  mildly 

accelerating  flow  over  the  plate)  and  the  source-rate  densities  which 

are  also  dependent  on  U . Once  the  form  of  the  source- rate  density, 

g,  is  accurately  determined  from  controlled  laboratory  experiments,  the 

use  of  the  technique  utilizing  fixed  transducers  for  determining  the 

width  of  the  transition  region  will  be  possible  for  more  complicated 

situations.  For  example,  in  transition  experiments  on  moiels  or  full- 

scale  vehicles  where  the  streamwise  traverse  of  sensors  or  maintenance 
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of  constant  flow  ve  ocity  Is  difficult  and  the  number  of  fixed  sensors 
Is  limited,  the  technique  adopted  In  the  (resent  experiment  wll1  be 
useful  in  determining  the  transition  width. 

It  Is  concluded  that  .the  sou-ce  distribution  of  turbulent  burst: 
is  still  not  firmly  established.  The  present  results  tend  to  sub- 
stantiate Emmons'  model  for  predicting  burst  statistics  In  the 
transition  boundary  layer,  and  support  the  line  source  model  for  burst 
generation.  The  study  indicates  that  more  detailed  wind  tunnel  flat 
plate  measurements  are  still  needed  in  the  form  of  streamwise  surveys 
of  burst  statistics  at  constant  flow  velocities.  New  experimental 
data  obtained  throughout  the  entire  streamwise  extent  of  the 
transition  region  at  a fixed  flow  velocity  would  overcome  the  uncer- 
tainties of  interpretation  inherent  in  the  present  data,  which  was 
obtained  by  varying  flow  velocity  at  a fixed  streamwise  position. 
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PLAN  VIEW 


i 


Figure  2 Measured  Properties  Of  Wall  Pressure  Bursts 
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Figure  3 Intermittency  Distribution  Through  The  Transition  Region 


*Run  No.  relates  data  to  additional  flow  parameters  giver,  in  Referenc 


Figure  5 Average  Burst  Length  Versus  In termi t tency  Factor 


THEORY 


Figure  6 Comparison  Of  Experimental  Data  With  Emmons 


Table  1 Emmons'  Probabilistic  Model 
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APPENDIX  A 


DERIVATION  OF  THE  INCREMENTAL  VOLUME  ELEMENT  j 

3 | 

Schubauer  and  K!ebanoffJ  found  the  shape  of  turbulent  bursts  to 

s 

be  wedge-like  pointing  downstream.  The  shape  of  the  Retrograde  Cone 
shown  in  Figure  A1  also  forms  a wedge  like  area  when  projected  on  the 

! 

xy  plane.  The  locus  of  all  points  PQ  which  can  influence  the  state 
of  turbulence  at  point  P forms  the  volume  represented  by  the  Retro- 

» 

grade  Cone  in  Figure  Al.  If  the  x = constant  cross  section  of  the  cone 

I 

is  approximated  by  the  shape  of  a triangle  the  volume  element  dVQ  in 
the  integral  relations  (4)  and  (8)  can  easily  be  found. 

The  volume  of  the  Retrograde  Cone  for  bursts  originating  at  xQ 
is  given  as 

v-K 5 


where  E,  is  the  downstream  distance  from  the  measurement  point,  P,  to 
the  point  of  burst  origin,  P . The  area  Afa  is  the  projection  of  the 
Retrograde  Cone  on  the  yt  plane  at  point  PQ.  From  Figure  Al  it  is 
seen  that 

\ ' 1 tb  d 

max 


where  d = 2£  tan  ( «)  , thus  Ab  becomes 


^ - £xh  tan  (a) 
max 

The  relationship  for  can  be  shown  to  be 

bmax 


B u u 

max  c t 
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thus  ^ is  then  given  as 


u 


Afc  " S u u tan  (a) 


c t 

and  the  expression  for  the  volume  is 

:»  U 


1 

3 


y.±i  r .X  — El 
* V 3 U.U 


tan  (a) 


c t 


Now  defining  a propagation  parameter,  o*,  as 

u u 


a*  - uu  tan(ct) 

c t 


we  can  then  express  V as 


« £i  si. 


3 LL 


From  this  relationship  the  incremental  element,  dV,  of  the 
Retrograde  Cor.e,  can  be  expressed  as 

dv-  £ 
e 

In  this  relationship  the  £ variable  represents  the  distance 
downstream  from  the  transition  point,  thus  in  keeping  with  the  nomencla- 
ture used  in  the  text  the  variable  £ is  expressed  as  x-xQ,  thus  we 
get  that 


dv  - (x-x  )z  dx 
o U o o 
e 
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APPENDIX  B 


MEASUREMENT  OF  INTERMITTENT  PRESSURE  FIELD  PROPERTIES 

The  fluctuating  pressures  on  the  plate  surface  were  measured  with 
Bruel  and  Kjaer  1/4-inch  and  1/8-inch  Condenser  Microphones.  To 
facilitate  the  measurement  of  the  pressure  on  the  plate  surface  over  as 
small  an  area  as  possible,  the  microphones  we»e  used  with  solid  pro- 
tective caps  mounted  flush  with  the  plate  surface  and  in  each  of  which 
was  drilled  a single  1 /32-inch  diameter  hole. 

The  intermittency  factor  of  the  wall  pressure  field  was  varied  by 
changing  the  flow  velocity,  since  the  position  of  the  pressure  trans- 
ducers were  fixed  and  could  not  be  varied  without  stopping  the  wind 
tunnel.  It  is  assumed  that  properties  of  the  intermittent  pressure 
field  which  are  dependent  on  y,  are  independent  of  whether  the  intermit- 
tency was  varied  through  change  of  flow  velocity  U^  or  change  in  the 
streamwise  position  x in  the  intermittent  region.  Inherent  in  this 
assumption  is  the  basic  assumption  that  secondary  effects  on  the 
transition  process,  such  as  turbulence  intensity,  plate  vibrat*  t,  etc., 
are  not  altered  by  the  small  fractional  change  in  velocity  necessary 
to  vary  the  intermittency  factor  over  its  range  from  0 to  1 , at  a fixed 
point  on  the  surface.  The  fact  that  the  distributions  of  the  inter- 
mittency factor  through  the  intermittent  streamwise  transition  region 
collapsed  well  when  represented  in  terms  of  a normalized  Reynolds 
number  (as  shown  in  Figure  3)  further  substantiates  the  interchange- 
ability  of  velocity  or  streamwise  position  for  this  study  of  the 
intermittent  pressure  field. 
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The  intermittent  signals  from  the  pressure  transducer  located  in 
the  intermittent  transition  region  of  the  flat  plate  boundary  layer 
were  displayed  on  oscillographs.  By  noting  the  fraction  of  the  sample 
time  that  the  oscillograph  trace  was  turbulent  the  intermittency  factor 
was  determined  as  a function  of  flow  velocity  and  streamwlse  distance  x 
from  the  plate's  leading  edge.  The  average  burst  length  and  burst 
frequency  were  also  calculated  from  the  oscillographs.  Simultaneous 
oscillograph  traces  of  the  signals  from  two  transducers  located  on  the 
plate  surface  with  streamwlse  spatial  separations  were  used  to  deter- 
mine the  velocities  of  the  leading  (downstream)  and  trailing  (upstream) 
edges  of  the  pressure  bursts.  Simultaneous  oscillographs  of  signals 
from  transducers  separated  in  the  transverse  direction  to  the  mean 
flow  were  used  to  determine  the  width  of  the  bursts  and  the  interior 
wedge  angle  0 as  shown  in  Figure  2.  These  properties  of  the  pressure 
bursts  determined  from  simultaneous  signals  of  two  spatially  separated 
microphones  were  accomplished  by  measuring  the  time  differences 
between  the  detections  of  the  bursts  at  the  two  transducers. 


APPENDIX  C 


METHOD  OF  NORMALIZATION  OF  THE  DATA 
CASE  1 - DIRAC'S  DELTA  FUNCTION 

The  statistical  relationships  which  have  been  derived  using  a Dirac's 
Delta  Function  as  the  assumed  form  of  the  source  rate  density  function 
are 
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The  streamwise  coordinate  used  in  these  equations  can  be  selected  such 
that  transition  occurs  at  x = 0,  thus  xt  - 0 and  = x.  This  allows 
the  streamwise  coordinate  to  be  normalized  over  the  transition  region 
as 


Xi  ” (Ax  ) 
t i 


(Cl) 


Equation  (21)  gives 


(Ax.) 
t i 
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thus 


xi " T^qjn 

y no* 
and 


0 < x < 1 


Using  Equation  (C2)  we  can  write 


x 


(2.046) 


and  Equations  (16)  and  (17)  become 


(C2) 


(C3) 


(C4) 


_a 

Yj  * 1 - e*P [-4. 185  xj 


£b  fin  J 

(A^r‘  nt«i(a)  2.046  exp[-4.185  I ) 


(C5) 

(C6) 


Equation  (C5)  Is  now  In  normalized  form.  Table  C1.0  shows  the  values 
of  y calculated  using  Equation  (C5)  for  values  of  0 _<  x1  <1.1. 

In  order  to  normalize  Equation  (C6)  we  define  a normalized  burst  rate 
frequency  fg  as 


_ 2 

2,046  x^  exp[-4.185  x^ ] 


(C7) 


Once_£Jie  relatloishlp  between  n and  Ug  is  known  fg^  can  be  divided  by 
n tan(=)  to  normalize  the  experimental  data  for  comparison  with 
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the  analytical  results  given  by 


f_  - 2.046  x exp [-4. 185  x2] 
0,  I 1 


Table  C1.0  shows  the  analytical  results  found  using  Equation  (C8). 
A relationship  between  n and  Ue  is  determined  by  setting 


nl/^*  tan(0° 

equal  to  the  maximum  burst  rate  frequency  calculated  from  Equation 
(C8).  Hence 


Then  we  can  write 


■ JS-- 


tan (a) 


■ .430  r 


tan (a) 


[0,430  tan(a)]; 


Using  the  experimentally  determine  values  of  tan  (<*)  and  a*, 

tan  ( <*)  = 0.571,  a*  = 3.327  tan  ( «)  = 1.898,  Equation  (C12)  becomes 


n ■ 31.165  ^ax 
U 
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With  the  experimentally  measured  values  of  fg  at  each  flow  velocity, 

Dmax 

U , it  is  then  possible  to  determine  the  relationship  between  n and  U . 
e c 

This  relation  between  n and  Ug  Is  shown  graphically  in  Figure  8.  For 
simplification  a linear  relationship  between  n and  Ua  Is  assumed  (as 
shown  by  the  straight  line  In  Figure  8)  for  normalizing  the  experimental 
values  of  fg  at  a given  Ug. 

Using  the  relationship  between  n and  U the  experimental  values 

/ 7C 

of  fg  are  normalized  In  the  form  (tan(«))  and  collared  to 

the  normalized  analytical  results  found  from  Equation  (C8).  This  com- 
partlson  Is  shown  In  Figure  6. 

CASE  2 - CONSTANT  SOURCE 

The  normalization  process  for  the  Equations  where  the  source  rate 
density  function  Is  a constant  Is  the  same  as  In  Case  1.  The 
statistical  relationships  are 


Y2  - 1 - £*] 

e 

fB  - fB  “ 8t*n(®)  exp^3u^"  { * C exP^ «3C 
2 e >q  e 

T - li*l  < T - 1*^? 
tb  fB(x>  - tb2  f^Tx) 


(23) 

(24) 

(25) 


The  normalization  process  of  the  streamwlse  coordinate,  x,  Is  similar 
to  that  In  Case  1.  In  this  manner  the  relation 


(Ax  ) - 1.448 
t 2 


(27) 
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is  used  to  yield  a relation  for  the  normalized  coord 1 * x . 

7 2 9 


(C14) 


where  again 


0 < x < 1. 
~ 2 “* 


(C15) 


Using  the  normalized  coordinate  In  Equations  (23)  and  (24)  gives 


Y2  - 1 - exp [-3. 035  xj  (C16) 

U 2/3  (x  , 

fB  “ 4‘361  a (ten  (a))  [^]  I 2 ; exp[3.035  U3-x  )]  d;  (C17) 

2 • 2 i 

O 

Equation  (C16)  is  now  In  normalized  form.  Table  Cl.l  shows  the  values 
of  y calculated  using  Equation  (C16)  for  values  of  0 < x < l.l. 

To  normalize  Equation  (C17)  we  define  a normalized  burst  rate  frequency 


fB  ' 


*“<«>[**  ] 


U_  2/3  “ 4,331  j 


fX2 


S exp[3.035(£3-x  )]  d? 


(C18) 


The  right  hand  side  of  Equation  (C18) 


- 4.361  J 


2 C exp(3.035(?3-x3)]  dC 

2 


(C19) 
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is  the  normalized  analytical  result.  Table  Cl.l  shows  the  analytical 
results  found  using  Equation  (Cl 9).  The  left  hand  side  of  Equation 
(Cl 8)  gives  a relation  for  normalizing  the  experimental  results. 


u 273 

gtan  (a) 


(C20) 


Once  the  relationship  between  g and  Ug  is  known.  Equation  (C20)  can 
be  used  to  normalize  the  experimental  data  thus  allowing  a direct 
comparison  between  the  normalized  analytical  results  (Equation  (Cl 9) ) 
and  the  normalized  experimental  results  (Equation  (C20) ) . 

A relationship  between  g and  Ug  is  obtained  by  setting  Equation 
(C20)  equal  to  the  maximum  burst  rate  frequency  calculated  from 
Equation  (C  19) 

f 


max 

0 TPT 


0.618 


(C21 ) 


gtan  (a)  [^*] 

From  Equation  (C21)  we  can  write 


£• 

g - 2max 

U2 (0.618  tan (a) ) 1 
€ 


(C22) 


then  using  the  experimentally  determined  values  of  tan  («)  and  o*  we 
get 


g 


82.202 


B2 


max 


(C23) 
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With  the  experimentally  measured  values  of  fR  It  Is  therefore  pos- 

max 

slble  to  determine  a relationship  between  g and  Uft.  This  relationship 
between  n and  Is  shown  In  Figure  9. 

Knowing  the  relationship  between  g and  Ufi,  It  Is  then  possible  to 
normalize  the  experimental  data,  using  Equation  (C20),  and  compare  It 
to  the  normalized  analytical  results  found  from  Equation  (C19).  This 
comparison  Is  shown  In  Figure  7.  Again,  a simple  linear  relationship 
fchown  on  the  figure)  was  used  for  normalizing  the  experimental  values 
of  fg  at  a given  Ufi. 
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TABLE  C1.0 


| 

f 

I 

i 

t 

i 

i 

1 


.05 

.10 

.15 

-20 

.25 

.30 

.35 

.40 

.45 

.50 

.55 

.60 

.65 

.70 

.75 

.80 

.85 

.90 

.95 

1.00 

1.05 

1.10 


NORMALIZED  ANALYTIC  RESULTS 
g * n6(x-xt) 


Yi 

f8i 

.01041 

.10122 

.10282 

.04986 

.19619 

.20892 

.08987 

.27928 

.32178 

.15414 

.34608 

.44539 

.23015 

.39372 

.58456 

.31384 

.42110 

.74529 

.40110 

.42881 

.93538 

.48809 

.41889 

1.16520 

.57150 

.39446 

1.44880 

.64875 

.35928 

1.80570 

.71803 

.31725 

2.26329 

.77834 

.27207 

2.86077 

.82935 

.22691 

3.65500 

.87135 

.18423 

4.72975 

.90502 

.14573 

6.21027 

.93133 

.11239 

8.2868 

.95128 

.08455 

11.25257 

.96629 

.06207 

15.56734 

.97711 

.04449 

21.96296 

.98478 

.03114 

31.62414 

.99009 

.02129 

46.50084 

.99368 
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.01422 

69.85757 

TABLE  Cl .1 


NORMALIZED  ANALYTIC  RESULTS 
g a constant 


A 

2 

Y2 

fB2 

tb2 

.05 

.00036 

.00545 

.06961 

.10 

.00303 

.02176 

.13925 

.15 

.01019 

.04876 

.20904 

.20 

.02399 

.08596 

.27910 

.25 

.04632 

.13247 

.34967 

.30 

.07869 

.18689 

.42103 

.35 

.12203 

.24722 

.49361 

.40 

.17656 

.31101 

.56771 

45 

.24164 

.37531 

.64385 

.50 

.31575 

.43697 

.72258 

.55 

.39650 

.49286 

,80449 

.60 

.48089 

. 54020 

.89022 

.65 

.56552 

.57682 

.98040 

.70 

.64695 

.60145 

1.07565 

.75 

.72212 

.61378 

1.17652 

.80 

.78863 

.61450 

1.28337 

.85 

.84500 

.60513 

1.39634 

.90 

.89061 

.58779 

1.51518 

.95 

.92591 

. 56484 

1.63925 

1.00 

.95194 

.53859 

1.76748 

1.05 

.97022 

.51106 

1.89844 

1.10 

. 98240 

.48382 

2.03052 
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